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Abstract 

The occurrence of mesoscopic fluctuations in statistical systems implies, 
from the point of view of dynamical theory, the existence of local instabil- 
ities. However, the presence of such fluctuations can make a system, as a 
whole, more stable from the thermodynamic point of view. Thus, in many 
cases, a local dynamic instability is a requisite for the global thermody- 
namic stability. This idea is illustrated by several spin models. 
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1 Introduction 



There are two general types of fluctuations in statistical systems, microscopic 
and mesoscopic. The former are small oscillations about a ground state, and 
they define collective excitations, such as phonons, magnons, etc. This kind of 
fluctuations is pertinent to equilibrium state. Collective excitations characterize 
a set of quantum states of a statistical system, and the word microscopic reflects 
the microscopic nature of such fluctuations. 

Contrary to these, mesoscopic fluctuations are such that make a macroscopic 
system heterogeneous or locally drive a system out of equilibrium. Locally can 
mean in space, or in time, or both. So, mesoscopic fluctuations are, in general, 
nonequilibrium and are related to an averaged description of a statistical system. 
The word mesoscopic, in relation to space, means that the characteristic size of 
such a fluctuation is much larger than the average interparticle distance but is 
much smaller than the size of a system itself. In relation to time, mesoscopic 
implies that the characteristic lifetime of a mesoscopic fluctuation is much longer 
than an effective oscillation period of microscopic fluctuations but much shorter 
than the observation time. When mesoscopic fluctuations correspond to the for- 
mation of nuclei of one thermodynamic phase inside another, they are called 
heterophase fluctuations [1]. The latter are ubiquitous in nature, and plenty of 
examples are described in review [1]. Recently, much attention has been paid to 
the study of mesoscopic fluctuations in high-temperature superconductors [2-7], 
where the corresponding phenomenon is often termed phase separation [3-6], al- 
though it would be more correct to call this effect mesoscopic phase separation 
in order to distinguish it from the principally different Gibbs phase separation 
occurring at macroscopic scales (see discussion in [1,7]). 

When mesoscopic fluctuations are frozen in time, the corresponding system 
looks like an ensemble of clusters with different properties. Such a system can be 
even equilibrium [8,9]. And if these fluctuations are not frozen in time, they make 
the system nonequilibrium [1]. The latter case is more interesting than that of a 
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frozen macroscopic structure, since it involves three difficult questions: (i) How to 
develop a statistical description of such a nonequilibrium and nonuniform system? 
(ii) Can this mesoscopic state be an attractor and, if so, what kind of attractor is 
it? (iii) Though the existence of nonequilibrium mesoscopic fluctuations implies 
local instability, but is it possible that globally the system is nevertheless stable? 

The answer to the first question has been done by developing a consistent 
statistical theory of systems with such mesoscopic fluctuations [1]. In this ap- 
proach, after averaging over these fluctuations, a renormalized Hamiltonian is 
defined representing a set of phase replicas, each of which describes an effective 
equilibrium system. To deal further with the renormalized Hamiltonian, one may 
employ the techniques of equilibrium statistical mechanics. 

As an answer to the second question, it has been conjectured that the state of 
a system with nonequilibrium mesoscopic fluctuations is a chaotic attractor [1] . It 
was also shown [10] that a uniform statistical system is structurally unstable with 
respect to arbitrary small random external perturbations. One can recollect as 
well that in dynamical theory there are plenty of examples of dynamical systems 
with chaotic attractors (see e.g. [11-13] and references therein). 

The third question has been considered for some simple spin models [1] and 
analyzed in more detail for a chaotic lattice-gas model [14]. In the present paper, 
the study of several other less trivial spin models is given, with the aim to show 
that for each model there can be found a region of parameters, where mesoscopic 
fluctuations do make the system globally stable, that is, thermodynamically more 
stable than the analogous system without these fluctuations. 

2 Effective Hamiltonian 

This section is a very brief recollection of the main definitions we shall need in 
what follows for analyzing systems with mesoscopic or heterophase fluctuations. 
We employ the theory of such systems developed in Ref. [1]. Throughout the 
text, the terms mesoscopic and heterophase will be used in parallel. 
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Imagine that we are dealing with a system in which randomly in space and 
time there arise mesoscopic fluctuations. For crystals and liquids, these could be 
fluctuations of local space structure [15-18]. In the case of spin systems, these 
are fluctuations of local magnetization [1]. Following the general theory [1], we 
can average over such stochastic mesoscopic fluctuations and obtain an averaged 
effective Hamiltonian 

H eff = @ V H V , (1) 

in which the index v enumerates qualitatively different phases, and H u is a phase- 
replica Hamiltonian representing a pure //-phase. Each Hamiltonian H u is defined 
on a Hilbert space Ti u of microscopic states typical of the corresponding //-phase 
[1,19,20]. The averaged Hamiltonian (1) acts on a fiber space 

y = ®uH u . (2) 

Hamiltonian (1) depends on a set {w u } of geometric phase probabilities w u , with 
the properties 

X>„=1, 0<w v <l. (3) 

These phase probabilities are defined as quantities providing an absolute mini- 
mum for the free energy 

/ = -llnTrexp(-/?tf e// ), (4) 

where T is temperature, N is the averaged number of particles, (5T =1, k B = 1, 
and the trace is taken over the fiber space (2). The phase probabilities as functions 
of thermodynamic and Hamiltonian parameters are given by the solutions of the 
equations 

df d 2 f 

0, ^>0, (5) 



dw v dwl 
under the normalization condition (3). 

As concrete examples, we analyse below spin systems. Heterophase states 
in these systems are characterized by mesoscopic fluctuations of local magne- 
tization. Such spatial fluctuations can be observed in experiment by means of 



neutron scattering methods, like diffuse scattering, small-angle scattering, Bragg 
reffections, and polarized-beam scattering [21]. 



3 Model with Competing Interactions 

Consider a one-dimensional Ising-type model modified by including two kinds of 
spin interactions, long-range and short-range interactions, 

Jij = alS^j^ + (1 - a) Jg, (6) 

where J°- satisfies the properties 

and a is a crossover parameter [22,23]. For a = 0, we have only a long-range 
interaction, while for a — 1, we get the Ising nearest-neighbour interactions. 
Assume that mesoscopic fluctuations are heterophase fluctuations between two 
phases, ferromagnetic and paramagnetic. The averaged Hamiltonian (1) consists 
of the phase-replica terms 

Hu = w 2 u {^NU - 1 £ JijSis)j , (7) 

where iV is the number of lattice cites, U is a crystalline-field parameter, and 
Si = ±1. Note that the crystalline field cannot be omitted since it influences 
the values of the phase probabilities, although it does not contain spin variables 
[24] . This makes the situation rather different from the case of pure monophase 
systems. 

Each phase is characterized by an order parameter 

1 N 

a ^= J?Y,( S i)^ ( 8 ) 
i=l 

where (. . ) v implies the statistical averaging over the space of typical states H u ; 
that is, an average {A) u of an operator A means 

exp(-(3H u ) 



{A) u = Tr v p v A, p v = 



Tr„exp(-PH V ) ' 



where Tr u means a trace over Ti v . Let v — 1 correspond to ferromagnetic phase 
and v — 2, to paramagnetic phase. Then, by definition, 

a x £ 0, a 2 = 0. (9) 

This condition relates the order parameters (8) with the space 7i v of typical 
states, making it possible to construct the latter as quantum weighted spaces [1]. 

Calculating the specific free energy (4), we can use for the short-range part of 
the Hamiltonian the transfer-matrix method (see e.g. [25]), and the long-range 
part, as is known [26], is asymptotically equivalent to the mean-field form. As a 
result, we obtain 

1\ „ w 2 



f = ( w 2 - w + Ij U - ^- aT-{\- a)Ja< 



-Tin 



cosh(/9 + \J sinh 2 (/9 + exp(— 4</?i) 



Tln(2cosh^ 2 ) , (10) 



where u> = u>i, <r = <7i, and 



_ 2 (l-«)^o- _ 2 al _ 2 al 

¥ = W — , <p 1 = W—, <p 2 = (l-w)—. 

For what follows, it is convenient to introduce dimensionless quantities 

U I T 

u =r 9 =r 

For the probability of the ferromagnetic phase, from the first of equations (5), we 
find the equation 

4gwaexp(— 4</?i) 
coshy^sinh 2 ^ + exp(— 4<pi) + smh 2 ip + exp(—Aipi) 

+ u(2w - 1) - 2agw - 2(1 - a)wa 2 + 2a(l - w)gtanhcp 2 = 0. (12) 

For the order parameter a = oi, defined in (8), we get 

sinho? 

/ 2 = • (13) 

y'sinh ip + exp(— 4(/?i) 

The heterophase state is thermodynamically more stable than the pure state, if 
the inequality 

A/ = /(l) -/(«;)> (14) 



holds true. When (5) is valid but (14) is not, the heterophase state is metastable. 

Let us investigate the stability of the system at T = 0. Then, from (12), we 
obtain the ferromagnetic-phase probability w(T = 0) = Wq in the form 



Also, we have 



2u-\a\g 

w o — - A ; 7 i — — • (15) 

Au - 1 + a - (a + a \)g 



d 2 f 

2(4w — ag — 1 + a — \a\g), 



dw 2 

Af _ (2it - ag - 1 + a) 2 (4u + 2ag - 1 + a 
~T ~ 4(4m - 1 + a) 2 

Af (2u-ag-l + a) 2 



a < 0, 



a > 0. 



J 4(4u-2ag- 1 + a)' 
From conditions (3), (5), and (14), it follows that the heterophase state is 
absolutely stable if either 



u 



> max — a — 2ag), -(1 — a + ag)^ , a < 0, (16) 



or 

u>hl-a + ag), a>0. (17) 
The heterophase state is metastable if either 

-(1 - a + ag) < u < ^(1 - a - 2ag), a < 0, (18) 



or 



In the case when 



u < -ag, a > 0. (19) 



^\a\g < u < -(1 - a + ag), (20) 

the system at zero temperature is purely ferromagnetic, but beginning from a 
finite temperature T n , called the nucleation temperature [1] and defined by the 
condition w(T n ) = 1, the heterophase state becomes profitable, being a mixture 
of ferromagnetic and paramagnetic phases. 
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At low temperatures, when t — > 0, we find the following asymptotic behaviour 
for the ferrophase probability 

w 4(1- a + ag) f i»L , x l 

— ~ 1 - - yi , , exp <^ — 1 - a + Aag) \ + 

wq Au — ag — 1 + a — \a\g [ t J 

2\a\g(2u - ag - 1 + a) j (1 - w ) 2 \ 

+ — ; — p— 7- ; — exp < 2\a\g } , 

{2u — \a\g)(Au — ag — 1 + a — \a\g) { t J 

the order parameter 



a ~ 1 - 2 exp < — ^(1 - a + Aag) > , 



the entropy 



S ~ -a + «5f)exp| — ^(1 - a + Aag) | + 

, (l-^o) 2 , , / (l-^o) 2 ol , \ 
+ 2t l«l^ ex P{ — 2 I^)> 

and for the heat capacity 

C y ~ -^(1 - a + a#) 2 exp ^ — ^(1 - a + Aag) \ + 
(1 - w ) 4 2 2 f (1 - W 2 



a g exp < 2\a\g 



■ ^2 - - ~^ ^ t 

The positivity of the specific heat indicates that the heterophase state is stable 
with respect to thermal fluctuations. 

Now, let us analyse the critical behaviour of the model. The critical temper- 
ature t c is defined by the condition a(t c ) = 0, which gives 

1 — a / aq\ 
^ = —5— exp hf • (21) 



As follows from (21), there exists a negative value of the crossover parameter 
a = cko, 

a = l — {eg > 1), (22) 

eg - 1 

such that for a < ao the ferromagnetic state is impossible at all temperatures. 
But if g < e _1 = 0.3679, then a positive solution for t c is available for any 



a < 1. The appearance of the limiting value (22) is quite explicable. Really, 
negative values of a correspond to the antiferromagnetic character of the short- 
range interaction. The presence of an interaction having the opposite sign, as 
compared to the ferromagnetic long-range interaction, serves as a disordering 
factor. The onset of ferromagnetic order is possible only it the disordering short- 
range interaction is not too large. One might recollect several other examples 
when an ordering in a system occurs only if some limiting relations between 
competing interactions take place. Recall, for instance, the criteria of magnetism 
in the Hubbard model [27,28]. 

It is interesting that the crossover behaviour of the critical temperature, be- 
tween the mean-field value t c — | (a — 0) and the short-range case t c — (a — 
1), is nonmonotonic. The maximum of (21) occurs at 

, _ 9 n ln ff / 9 ox 

tmax - B ( 1 + hlg y a max 1 + ]ng ' W 

The ratio of this maximum to the mean-field critical temperature t c — 1/8, that 
is, 

Rt 9 



l + m#' 

can become arbitrary large for g ^> 1. 

The presence of mesoscopic fluctuations, as well as the antiferromagnetic 
short-range interaction, can change the second-order phase transition to the 
first-order one. The region of first-order phase transition is defined either by the 
inequalities 

u < u < u t (a > (Xi), (24) 

or by the inequalities 



where 



u < u < u — \u t — u \ (a < a < oil), (25) 
12(1 -«)(«£■ + 2£ c ) 2 



u t = u + 



m - (i - «) 2 



Aa{l-a)g + Atl-{l-af 
U ° = 4(l + t c -a)* ' (26) 
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e*i = 1 - 2VSt c . 

The value u = u t , denned in (26), is a function of the parameters a and g. The 
equation u = u t (a,g) describes a surface on which the order of phase transition 
changes. This is called a tricritical surface. On the latter, the critical indices also 
change by a jump. Thus, considering the asymptotic behaviour of the specific 
heat Cy oc \r\~ a , order parameter a oc \t\^ , and susceptibility x °^ M~ 7 ? when 
approaching the critical point t c , so that r = (t — t c )/t c — > —0, we obtain 

a = 0, p = ± 7 = 1 (u^u t ) (27) 

outside the tricritical surface, and 

<*=\, P = \, 7 = 1 (u = ut) (28) 

on the tricritical surface. 

One more critical index can be introduced [1] for the phase probability w as 

w-ioc|r| £ . (29) 

This index e is specific for heterophase systems. In our case we find that this 
index also jumps on the tricritical surface u = u t (a,g), 

1, u ^ u t 

' T (30) 

2, u = U t - 

Despite the seeming simplicity of the model considered, it displays quite non- 
trivial behaviour. One of the most interesting features is a strongly nonmonotonic 
dependence of the critical temperature on the crossover parameter a. Also, in 
the space of three parameters, u, a, and g, there is a region, where mesoscopic 
fluctuations make the system thermodynamically more stable than that system 
without such fluctuations. 



4 Heterophase Spin Glass 



Consider a generalization of the Sherrington-Kirkpatrick spin-glass model [29] 
to the case of a system with mesoscopic fluctuations. We shall keep in mind 
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paramagnetic fluctuations inside the spin-glass phase [30]. Let interactions J^ 
between spins at sites i,j — 1, 2, . . . , N be distributed by the Gaussian law 



1 fN\ 1/2 f N / J 



2~ 



F «'>=jfeJ eXP |-2^l J «-Fj }• (31) 

For simplicity, we put Jo = 0. The average over interactions for a function A{ J^} 
of a set { Jij} of interactions J^ is defined as 

[A{J tJ }] av = J A{J l3 } J] P( J«)^« • (32) 

For an Ising-like system with mesoscopic fluctuations, following the general 
renormalization procedure [1], we have an effective Hamiltonian (1) with the 
phase-replica terms 

HMi) = w 2 v {^NU - E Ji jSi s}j , (33) 

in which U is a crystalline-field constant and Sj = ±1. To distinguish phases, 
we need an order parameter. For spin glasses, this is the Edwards-Anderson [31] 
order parameter 

9, = , (34) 

where (. . .) v means a statistical averaging with the Hamiltonian (33) over the 
space 7i v of states typical of a phase with the parameter (34), and [. . ] av denotes 
the averaging (32) over interactions. Let the spin-glass phase be indexed by v — 1 
and the paramagnetic phase, by v = 2. Then, by definition, 

qi ^0, q 2 = 0. (35) 

The ferromagnetic phase, because of Jo = 0, is absent. 
The free energy (4) writes 

f = fi + / 2 , /, = [U{Jii}\ av , (36) 

where 

T 

UUij} = -—]nTr v exp(-PH v {J i: j}), 
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and Ti v = {si\ i = 1,2, . . . , N; q v } is the space of states typical of a phase with 
the order parameter (34). With the replica trick [31], one has 

[\nZ{J t A] = lim - ( [Z n {J i7 -i] - 1) = lim ^- [Z n {J l7 }] . 

Using this, for the free energy (36), we obtain 



U 



1 i 

" '2.. .4/i „\2 1 n T 2/i „..\4 



/ = - w 2 + (1 - w) 2 - -(3JW(1 - qY - -3J\1 - w) 



2 L J 4 4' 

/+oo . , 

p(x) In 2cosh(/3 Jw 2 q 1/2 x) dx-T\n2, (37) 
-oo L 

where w = w±, q = q\, and 

p{x) = vb exp (-^ 2 )- 

For the spin-glass order parameter (34), with v — 1, we get 

/+oo . . 

p(a;)tanh 2 (/3Jw 2 g 1/2 a;) tte. (38) 
-oo ^ ' 

The trivial solution q = of (38) is excluded by condition (35). From the first of 
eqs. (5), we get an equation 

w 3 (l-q) 2 -(l-w) 3 -u(2w-l)t = (39) 

for the probability w of the spin-glass phase, where 

U T 
- 7 , *- 7 . (40) 

Among three roots of eq. (39), we have to choose that one satisfying the normal- 
ization condition (3). In comparison to the standard spin glass [31,32], the phase 
probability w plays the role of an additional order parameter. 

Let us analyse the thermodynamic characteristics of the heterophase spin 
glass. It is convenient to introduce the notation 

u = 2W^ = 1.595769. (41) 
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At low temperatures, the spin-glass order parameter behaves as 
q - 1 - y£ - Uo{u - u)t 4/3 (u<u ), 

q^l->t--t 2 («>«o), (42) 

Z 71 

where t — > 0. The asymptotic, as t — > 0, behaviour of the phase probability is 

W ~ 1 — (m — ■u)t 1 / 3 (-U < Wo)) 

w ~ 1 (u > u ). (43) 
For the specific heat and entropy, we find 

Cy ^ 1 -{U -Uflh- 2 I* (U<U ), 

Cy ~ 1 - W (u > u ), (44) 

Z47T 

as t — > 0, and, respectively, 

5^-i(«-«o) 4/3 r 2/3 («<«o), 

5 ~ In 2 - — = 0.53399 (u > u ). (45) 

27T 

As follows from these expressions, the ground state is a pure spin-glass phase: 
w — 1 at t — 0. 

When the crystal-field parameter tt < u , the system is unstable at low tem- 
peratures, since Cy — > oo and 5 — > — oo, in analogy with the Sherrington- 
Kirkpartick case [29]. But for u > Uq, mesoscopic paramagnetic fluctuations 
stabilize the system making the behaviour of the specific heat and entropy nor- 
mal. 

In the vicinity of the critical point 

tc = \ (? = 0), (46) 
we have for the spin-glass parameter 

t-t c 

c 



q~\r\ [r = —^^-0) (47) 
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and for the phase probability 

1 r 2 

w~-- - . (48) 

2 4(u-3) y J 

From the stability condition 

^-L~2J( M -3) >0 (*-f c ), (49) 

we conclude that the second-order transition between the spin-glass and param- 
agnetic phases occurs if u > 3. The value u — u t — 3 corresponds to a tricritical 
point, where the second-order transition changes for the first order transition 
being realized for u < 3. The critical index for the phase probability, defined in 
(29), is 

In w — \ 

e = lim ; , , =2 (u < 3). (50) 



t-*o In \t\ 

In this way, mesoscopic paramagnetic fluctuations, when u > uq, stabilize the 
Sherrington-Kirkpatrick mean-field spin glass, making its specific heat finite and 
entropy positive. In order to check whether the heterophase spin glass becomes 
absolutely stable, one has to consider as well the sign of magnetic susceptibility. 
According to our analysis [30], the latter is positive, at least in the critical re- 
gion. The mean-field glass, as is known, can be made stable by invoking, for the 
Edwards-Anderson order parameter, solutions with a broken replica symmetry 
[33]. Generalizing this type of spin glass to the case including mesoscopic fluc- 
tuations, it is possible to show [30] that, again, for sufficiently large crystal-field 
parameter u, the free energy of the heterophase spin glass becomes lower than 
of a pure spin glass. That is, mesoscopic fluctuations can make the spin-glass 
system thermodynamically more stable. 



5 Systems with Magnetic Reorientations 

The appearance of coexisting magnetic phases with different directions of mag- 
netization is characteristic of spin-reorientational transitions in small or zero 
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external magnetic fields. The existence of such mixed states is well documented 
by a large number of experiments and have been discussed in detail in books 
[34,35] and reviews [1,36,37]. The standard theoretical description of magnetic 
reorientations, going through intermediate mixed states, is done by means of Lan- 
dau expansions involving a set of fitting functions taken from experimental data 
[34,35]. This is a purely phenomenological treatment giving no physical insight. 
Another approach to systems with magnetic reorientations can be based on the 
theory of mesoscopic fluctuations [1], with taking account of phase fluctuations 
corresponding to different angle phases with mutually orthogonal magnetizations 
[24,38,39]. Following such a microscopic approach, it is possible to show that the 
balance between phase probabilities and, as a result, magnetic reorientations are 
governed by the tendency of a system to reach the state of an absolute thermo- 
dynamic stability by allowing the appearance of mesoscopic fluctuations. 

Consider a system in which there can coexist four different phases, three of 
them, magnetic, having orthogonal to each other nonzero magnetizations and 
one, paramagnetic, with zero magnetization. These phases will be enumerated 
by the index v — 1,2, 3, 4, and the notation {u} = {a, 4}, where a = 1,2, 3, will 
be used. To separate the phases, we need an order parameter, whose definition, 
as usual for spin systems, is based on the average spin operator 



in which a — 1, 2, 3 and Si is a spin operator at a lattice site i — 1, 2, . . . , N. The 
vector order parameter of a //-phase is defined as 




N 



E Si= {s a } , 



(51) 



i=l 



V v = (<?)„ = {<}, 



(52) 



so that 



d = S, 



>a/3Va , 




(53) 



where 




(54) 
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In another way, we could write 



Va= 0, 



(55) 



where r] a is given by (54) and e a is a unit vector along the a-axis. The order 
parameters (55) define the directions of magnetization for the related phases, 
three of which are magnetic and one is paramagnetic. 

Following the general theory [1] for a spin system with anisotropic interactions 
J°!j, after averaging over random phase configurations, we come to an effective 
Hamiltonian (1) with the phase-replica terms 



where K is a mean kinetic energy per site for electrons and ions. Recall [24] 
that the term NU is the total potential energy, not including spin operators, 
of electrons and ions in a crystalline lattice. Therefore, the parameter U can 
be called the crystal-field parameter, structural constant, configurational energy 
per site, or lattice energy per site. In general, the values K, U and can also 
depend on the type of a phase, but, for simplicity, we assume that they are the 
same for all thermodynamic phases. 

Minimizing the free energy (4), we get an equation 




(56) 




(57) 



for the phase probabilities w u , where 



(58) 



a =l 



From the second of eqs. (5), we have the stability condition 




(59) 



while the inequality < w u < 1 is valid when 



U > sup{B a }. 



(60) 



a 
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For ferromagnets, the interactions are positive, thence B v > 0, and for antifer- 
romagnets, J° a are negative, hence B v < 0. Therefore, heterophase fluctuations 
appear easier in antiferromagnets than in ferromagnets, as follows from (59) and 
(60). 

To proceed further, we need to invoke some approximation. In what follows, 
we use the mean-field decoupling 

which yields 

2 2 1 N 

B a = J a Si] a , B 4 = 0; Ja = j^^Jij, 
where S is a spin value and r\ a is given by (54). For S — 1/2, we have 

i] a = ^tanh(f3w 2 a J a i] a ), 774 = 0. (61) 
The free energy (4) becomes 

/ = E { W K U + J ^l) - Tln[2cosh(/3^j a?7Q )]} + ^ _ Tln2 . (62) 

a=l 

In this way, the thermodynamic behaviour of the system is defined by the set 
of seven coupled equations: four equations for the phase probabilities (57) and 
three nontrivial equations for the order parameter (61). Among all admissible 
solutions, one has to choose those satisfying all stability conditions and providing 
an absolute minimum of the free energy (62). In addition to the case when the 
phase probabilities are found from (57), we need to consider the cases when one 
or several phase probabilities are put zero [38]. This implies that we have to 
compare fifteen types of solutions: 

(1) wi^0, w 2 ^0, w 3 ^0, w A ^0; 

(2) Wl ee 0, w 2 ^ 0, w 3 ^ 0, w 4 £ 0; 

(3) wi ^0, w 2 = 0, w 3 ^ 0, w 4 ^ 0; 

(4) vox ^0, w 2 ^ 0, w 3 ee 0, w A ^ 0; 
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5) wi = 0, W2 = 0, u>3 ^ 0, u> 4 ^ 0; 

6) Wi = 0, u> 2 ^ 0, w 3 = 0, u> 4 ^ 0; 

7) wi ^ 0, u> 2 = 0, w 3 = 0, u> 4 ^ 0; 

8) Wi = 0, u> 2 = 0, w 3 = 0, u»4 = 1; 

9) ^ 0, u> 2 ^ 0, u> 3 ^ 0, w 4 = 0; 

10) wi = 0, w 2 ^ 0, u>3 ^ 0, w 4 = 0; 

11) wi ^ 0, u> 2 = 0, w 3 ^ 0, u> 4 = 0; 

12) ^ 0, w 2 ^ 0, w 3 = 0, w 4 = 0; 

13) wi = 1, u> 2 = 0, w 3 = 0, u> 4 = 0; 

14) wi = 0, u> 2 = 1, w 3 = 0, u> 4 = 0; 

15) wi = 0, u> 2 = 0, w 3 = 1, u> 4 = 0. 

Leaving aside thermodynamics of this model [38,39], we would like to concen- 
trate here on the fact that choosing the most stable solutions makes it possible 
to describe various reorientation transitions. 

Let us arrange the exchange integrals in the order of their magnitude so that 

< Ji < J 2 < J 3 . 
The reorientation temperatures are defined by the equations 

Va (T a ) = (a = 1,2, 3). (63) 
The largest reorientation temperature is the critical temperature 

T c = supT a (64) 

a 

for a ferromagnet-paramagnet phase transition. When this transition is of first 
oreder, the transition temperature will be denoted by T . The temperature at 
which a pure state transforms into a mixture is called [1] the nucleation temper- 
ature. The latter is defined as 



0, 



(65) 



18 



where u — 1,2, 3, 4. 

The arising sequence of phase transitions can be classified according to the 
value of the parameter U corresponding to the potential lattice energy per site. 
When U < 0, no heterophase states appear, and no reorientation transitions 
occur. The sole transition is the ferromagnet-paramagnet phase transition at 
T c = \ J$ to which the notation 

[0 1] <- 2, T c -> [0 0] 

can be ascribed, where the numbers inside the square brackets means the existence 
of a nonzero or zero magnetization along an axis a = 1, 2, 3, respectively, and the 
number between the arrows shows the phase-transition order at a temperature 
T 

When < U < U , where the value U depends on magnitudes of J a , then 
there appears the sequence of phase transitions 

[0 1]^l,T n ->[lll]<-l,To->[0 0], 

in which T n = 7\ — T 2 . Both the nucleation and ferromagnet-paramagnet tran- 
sitions are of first order. 

Increasing further the lattice-energy parameter, in an interval Uo < U < U\, 
we get the sequence 

[0 1] <- l,T n -> [1 1] 2,T 2 -> [1 1 1] <- 1,T ^[0 0] 

of phase transitions, with T n = T\. At the value U = U\, the nucleation temper- 
ature T n corresponds to a tri critical point. 
For Ui < U < U 2 , we have 

[0 1] 2, T n -> [0 1] <- 2, T x -> [1 1] <- 2, T 2 -> [1 1 1] 1, T -> [0 0]. 

The nucleation becomes a second-order transition. 

When U 2 < U < U 3 , then the sequence of phase transitions simplifies to just 
one transition 

[0 1] <- 1,T ^[0 0] 
19 



between ferromagnetic and paramagnetic phases. For U = Us, the transition 
temperature T corresponds to a tricritical point. 
Finally, for U 3 < U < oo, we have the sequence 

[1 1 1] <- 2, T x -> [0 1 1] <- 2, T 2 -> [0 1] <- 2, T c -> [0 0]. 

As is seen, the existence of heterophase mesoscopic fluctuations makes it possi- 
ble to get a rich variety of reorientation magnetic transitions which are impossible 
in the pure case. Thus, mesoscopic fluctuations, at the same time, provide global 
stability for a system and make its physics much richer. Such fluctuations can 
also lead to specific features of thermodynamic characteristics, related to pre- 
transitional phenomena [1]. Probably, these fluctuations could be responsible for 
the existence of two length scales at structural and magnetic phase transitions 
observed with high resolution X-ray scattering techniques, particularly using 
synchrotron sources [40]. 
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